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Abstract 

■ The spectrum of core excitations of the Abrikosov vortex pinned by a nanohole of the size of 

^ . the coherence length is considered. While the neutral zero energy Majorana core state remains 

intact due to its topological origin, the energy of charged excitations is enhanced significantly 
compared to that in the unpinned vortex. As a consequence of the pinning the minigap separating 
the Majorana state from the charged levels increases from /S? /Ep for {Ep is Fermi energy, A - the 
bulk p-wave superconducting gap) to a signicant fraction of A. Suppression of the thermodynamic 
and kinetic effects of the charged excitations allows to isolate the Majorana state so it can be used 
for quantum computation. It is proposed that thermal conductivity along the vortex cores is a 
sensitive method to demonstrate the minigap. We calculate, using Butticker - Landauer - Kopnin 
formula, the thermal conductance beyond linear response as function of the hole radius. 
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Spin-triplet p-wave superfluids, both neutral, such as liquid He^ (and recently gener- 
ated by the Feshbach resonance on Li^ and K'^^) and charged such as superconducting 
material Sr2Ru04 and possibly heavy fermion UPt^ have resulted in very rich physics 
The condensate is described by a generally tensorial complex order parameter A exhibiting 
great variety of the broken symmetries ground states. The broken symmetry and boundary 
conditions give rise to the continuous configuration of the order parameter as nontrivial 

A 

topological excitations [2!]. Especially interesting is the case of the so-called topological 
superconductors, characterized by presence of electron-hole symmetry and absence of both 
the time-reversal and spin-rotation symmetry. Realizations of topological p-wave superfiuids 
are chiral superconductors like Sr2Ru04, with order parameter of the ± ipy symmetry 
type [3i] and the ABM - phase 1] of superfiuid He^ and other fermionic cold atoms [41] and 
topological superconductor CuxBi2Se3 that produces an equivalent pseudospin system on 
its surface ^. 

Magnetic field in type II superconductors easily creates stable line - like topological 
defects, Abrikosov vorticesj^, see Fig. la. In the simplest vortex the phase of the order 
parameter rotates by 27r around the vortex and each vortex carries a unit of magnetic fiux 
$0- Quasiparticles near the vortex core "feel" the phase wind by creating a set of discrete 
low-energy Andreev bound state. For the s-wave superconductors when the vortices are 
unpinned (freely moving) these states were comprehensively studied theoretically including 
the excitations spectrum [3], density of states j^, their role in vortex viscosity [sl, contribution 
to the heat transport [ol and to the microwave absorption [l^. The low lying spectrum of 
quasiparticle and hole excitations is equidistant, Ei = lu, where angular momentum I takes 
half integer values. The "minigap" in the s-wave superfiuids is of order of w = A'^/Ep 
<< A, where A is the energy gap and Ep is the Fermi energy. 

Free vortices in the p-wave superconductors exhibit a remarkable topological feature of 
appearance of the zero energy mode in the vortex core Uj . The spectrum of the low energy 
excitations remains equidistant, Ei = {l — l)u, but now / is integer 12|. The zero mode 
represents a condensed matter analog of the Majorana fermion first noticed in elementary 



particle physics 



13|]. Its remarkable feature is linked to the fact that its creation operator is 



identical to its own annihilation, called the state. Its topological nature ensures robustness 
against perturbations from deformations of order parameters and nonmagnetic impurities. 
The states obey non-Abelian statistics with their pairs constituting a qubit. This might 
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offer a promising metliod of tlie fault-tolerant quantum computation The main issue 
is to isolate the states from those above the minigap fl^ . 

The / = 1 Majorana mode in the p- wave vortex states has been a popular topic of 



study over the last several years 
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18|. While the minigap in the s and rf-wave 



superconductors was detected by STM [2l|], in p-wave it has not been observed. The major 
reason for that is the small value of the minigap u in the core spectrum (just mK for 
Sr2RuO^. To address this problem a current trend was to propose increasingly sophisticated 
combinations of materials and geometries. For example one of the proposals |20| to expose 
the Majorana state is to induce the s-wave superconductivity by the proximity effect on the 
surface of a topological insulator. The minigaps of the resultant non-Abelian states can be 
orders of magnitude larger than in a bulk chiral p-wave superconductor. 

In this note we propose to solve the Majorana minigap problem within the original system, 
a Px + iPy bulk superconductor in magnetic field by pinning the vortices on artificially 
fabricated dielectric inclusions of the radius comparable to the coherence length f of the 
superconductor [22]. It was shown theoretically [24] that in the s-wave superconductors 
pinning by an inclusion of radius of just R = 0.2 — 0.5^ changes dramatically the subgap 
excitation spectrum: the minigap u ~ A'^/Ep, becomes of the order of A. In the present 
note we present the spectrum and wave functions of the core excitations in the chiral p- 
wave superconductor. The charged states for i? = 0.1 — 0.4,^ are significantly pushed up 
towards A, so that they therefore interfere less with Majorana state. To expose the modified 
charged spectrum we propose to measure heat transport along the vortex axis, see Fig. lb. 
The temperature gradient between the top side and the bottom side drives heat along the 
vortex axisjs]. The exponential temperature dependence of thermal conductivity is very 
sensitive to the minigap. The temperature when the material undergoes isolator - thermal 
conductor crossover rises from A/ 100 for unpinned vortices to A/ 10 for inclusion radius 
R = 0.4e. 

We start with the Bogoliubov - de Gennes (BdG) equations for the Px+ipy superconductor 
in the presence of a single pinned vortex. The vector potential A in polar coordinates, r, ip, 
has only an azimuthal component A^p (r) and in the London gauge consists of the singular 
part = hc/2er and the regular part of the vector potential that can be neglected for a 
type II superconductor [l2]. In the operator matrix form for a two component amplitude 
the BdG equations read: 
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where for anisotropic dispersion 
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^ = (^) e^'^ (^V. - V,) +\ [iiV, - Vy) s (r) e^^] | , 

with A being the "bulk gap" of order Tc. The equations, possess the electron- hole symmetry. 
The Ansatz 
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converts them (for any I there are radial excitation levels denoted by n) into a dimensionless 
form, 
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with dimensionless energy Sik^n = Ei^^n/ Here distances are in units of ^. We chose 
the order parameter as a discontinuous function vanishing inside the core r < R and s = 
tanh(r) for r > R. (The question of justification of using this form often used instead 
of the fully self consistent approach was extensively studied in literature on both s-wave 



and non- conventional pairing js], 12|). In the clean limit BCS (applicable to SrRu20^ 
^ = hk±/m±A, where 

kl/2m± = Ep/h^ - kl/2m,, (5) 
and for given there is just one dimensionless parameter 



7 = l/2k^^ = m^A/2h'ki. 



(6) 



The Ansatz, Eq.([3]) was chosen in such a way that the equations become real. In the 
presence of a hole of radius R we assume that the order parameter profile is still accurate 
for r > R. In a microscopic theory of the superconductor-insulator interface, (see 23[), the 
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order parameter rises abruptly from zero in dielectric, where amplitudes of normal excitations 
f = g = 0, to a. finite value inside the superconductor within an atomic distance a from 
the interface, namely with a slope oc 1/a. This means that the boundary condition on the 
amplitudes is consistent with zero order parameter at the boundary point r = R — a in the 
self consistency equation. The sample will be cylindric with radius L, see Fig. la. 

Qualitatively the spectrum of a single vortex in a hollow disk with internal and external 
radii R and L consists of two Majorana states, several pairs of Andreev bound subgap states 
both quasiparticles and holes (related to each other by the electron - hole symmetry) and a 
continuum of states above a threshold at Tg. The spectrum was calculated numerically by 
using NAG Fortran Library Routine Document F02EBF. It computes all the eigenvalues, 
and optionally all the eigenvectors, of a real general matrix, for various values of the inclusion 
radius = ~ 0.9^, external radius L ~ 10 — 20 and the universal parameter 7 (kz) = 
l/2k±C, ~ 10^^ — 3 ■ 10^^. For sufficiently small /, in addition to the continuum of states 
above the superconducting gap with n > 1, there are bound Andreev states that correspond 
to the lowest n = 0. For / = 1 and various inclusion radii R there is the Majorana state, 
for which / (r) = g (r) near the "internal" surface" extending to the distance ^ into the 
superconductor, see Fig. 2a where the radial density, 

p{r) = 2nr{\f{r)f + \g{r)f), (7) 

is given for sufficiently large external radius of cylinder, L = 20 and 7 = 0.03. As was noticed 



in ref. [17||, the function oscillates with period l/kj_^. There is also the second Majorana 
mode on the "external" surface for which / (r) = —g{r), given in Fig. 2b. The energy of 
the both states is exponentially small eik^ oc due to tunneling between themlOj. The 



surface state wave function practically does not depend on the inclusion radius for all values 
of R considered, although it does depend on L. When one considers a small cylinder of the 
width of several coherence length like the one shown in Fig. la, the two Majorana modes 
start to overlap. All the other / = 1 excitations, n > 1 , are above the threshold and will 
not be considered and the index omitted. 

At angular momenta / 7^ 1 topology does not protect the energy (more precisely its 
absolute value). The wave function of bound states, is pushed out of its position near the 
vortex core when the inclusion is absent. The state that was localized low energy at i? = 0, 
becomes delocalized and approaches the threshold at i? = 0.2^ and eventually a.t R = OA^ 
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merges into the threshold. The energies in units of A for wide range of angular momenta / 
for R = (the blue points), 0.1.^ (green), 0.2^ (red), 0.4^ (brown) are given in Fig. 3. The 



R = line is well approximated for |/| < 10 by the semiclassical linear formula while for 
higher momenta it approaches the threshold along a universal curve that is independent of 
the inclusion radius. When the hole is present, the dependence on / is no longer monotonic. 
It first rises to a maximum at Z = — 1 and subsequently has a local minimum with energy 
about that of the I = state. Consequently this becomes a new minigap that becomes of 
order of the bulk gap already at relatively small inclusion radius. Beyond this minimum the 
curve approaches the universal R = spectrum, so that wave functions are located far from 
the inclusion. Beyond radius R = 0.4^ the effect of the minigap enhancement is saturated. 
Only the value of 7 = 0.03 is shown in Fig. 3 although dependence on it is very weak for all 
7 << 1. The enhanced minigap allows to utilize the Majorana states. 

Thermal conductivity is an effective tool to demonstrate the minigap due to activated 
behavior of the electron contribution. To calculate the quasiparticle contribution to ther- 
mal conductivity along the vortex cores when the upper side of the vortex line is held at 
temperature Ti and the lower side at temperature T2, see Fig. lb, we use a general ballistic 
(width of the film smaller than mean free path) Kopnin-Landauer formulaj^. The heat 
current at temperature lower than the threshold to continuum of states is carried by the 
bound core states (except the Majorana). For a single vortex it consists of the contribution 
of quasiparticles and holes, 1 = 2 ^ /; (T2) — Ii (Ti), 



dkz 
27rh 



Ik^ 



Elk. 



1 + exp {Elk jry 



dkz 

where the energy depends on transferred momentum along the field kz via 7, see Eq.([6]). 
The maximal value of kz is k^^^ = \f2jm^Ep. Since Eik, is monotonic the integral can be 
transformed into 

27rni + exp {EjT) ^ 2^ ) ' ^ ^ 

Since the temperatures are below the threshold, Ep was replaced in the upper limit of the 
integral by 00. Here lower limit of integration is energy for kz = and 

n (x) = vrVe - xV2 + X log (1 + e^) + Lis (-e^) , (10) 
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where Li is the polylog function. For small temperature differences the linear response can 
be used, 



h dl t s—^ 



2n ( 



(11) 



t J 1 + exp (eio/t) 

where t = T/A. The values of energies of the core states, presented in Fig. 3, therefore 
allow to calculate the thermal conductivity. In Fig. 4 the heat conductance of a single vortex 
line is given as function of inverse temperature (in units of A^^). While for unpinned vortex 
(the blue line) the crossover from thermal insulator to conductor is not well defined in the 
relevant temperature range A/100 < T < A/4, for the radius of the inclusion R = 0.1^, 0.2^, 
0.4^ (brown, red and green respectively) one observes a well defined crossover temperature. 

To summarize, the wave functions and the spectrum of the core states well pinned vor- 
tices in the chiral p wave superconductor in magnetic field was studied, see Fig. 3. The 
pinning by a nanohole of the order of coherence ^ is assumed. It is shown that while the 
neutral Majorana mode is largely intact, all the other Andreev bound states have their 
energies significantly lifted. The minigap (the energy difference with the Majorana state) 
is consequently increased from its value for unpinned vortices, A'^/Ep, to a fraction of the 
bulk superconducting gap A. In superconductor Sr2RuOi while the minigap for unpinned 
vortices is just (using 3| A = 2K, Ep = 10'^ i^) u = AmK, it would become about 0.2K 
for a nanohole of radius of = 65nm [3] ) R = 0.4^ = 25nm. This would suppress the 
thermodynamic and transport effects of the charged excitations and allow to observe the 
physics of the Majorana state. 

We propose to demonstrate the enhanced minigap by measuring thermal conductivity 
along the vortex direction. The activated behavior of the electron contribution per vortex 
is given in Fig.4 and compared to the phonon contribution that dominates at temperatures 
below Tm = A/10. Magnetic field B » Hd creates = SB/^q vortices, see Fig. lb, 
over area S, so that heat conductivity is k = where is the sample width. For 

B = 0.5T (between Hd and Hc2 for in Sr2RuO^, with R = 25nm] Lz = 70nm at T = 0.2K 
one obtains conductivity k, = 0.05W/ Km. The phonon contribution in ballistic regime 
{Lz << Iph , where Iph is the phonon mean path) is I25I] n cc CvphLz- Taking the phonon 
heat capacity C in Debye approximation, one obtains for phonon conductance per vortex, 

tiKphl^ = (A/e^)' LznA%/B, (12) 



where ©i3,and ua are the Debye temperature, density of the atoms. This dependence for 
various vortex densities presented in Fig.4 (dashed hnes) demonstrates that for temperature 
above the minigap value, the quasiparticle heat conductance clearly dominates. In particular 

for parameters = 50nm,B = 0.5T, a = 0.5nm,T^ = 0.2( here a is the inter-atomic 
distance). 

It should be noted that vortices in such an experiment (for the field cooling protocol) 
will be trapped by the holes rather than by point defects. Those initially not pinned by the 
holes can be effectively pushed out by a small bias current. Our experimental proposal does 
not require an ideal hexagonal lattice of identical holes. Periodicity does not play a role and 
the array may even be random. 

Acknowledgements. B.Ya.S. and I.S. acknowledge support from the Israel Scientific Foun- 
dation. 
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Figure captions 



Fig. la. 

A single vortex in type II superconductor in magnetic field pinned on an insulator insertion 
of radius -R ~ ^ parallel to the field. The radius of the superconducting disk is L >> The 
geometry used to calculate the spectrum of the chiral p-wave core states. 

Fig. lb. 

Heat flow through a superconductor in magnetic field along the vortex cores. The tem- 
perature difference between the bottom (T2) and the top (Ti) contacts leads to energy flow 
carried by both the neutral (Majorana) and charged core states. Numerous vortices are 
pinned in inclusions and might be arranged as a lattice. 

Fig.2a. 

The radial density p(r), Eq.(I7j), of the Majorana core mode localized on the inclusion 
(green area) of radius R = 0.2^ as function of the distance from the vortex center. The 
value of the only parameter characterizing the system is 7 = 0.03. The energy of the state 
exponentially small as function of the sample radius L. The wave function is pushed out of 
the center by the inclusion, but otherwise remains intact compared to the unpinned vortex. 
It extends several coherence lengths inside the superconductor. 

Fig.2b. 

The radial density p (r) as function of the distance from the vortex center, Eq.Q, of the 
Majorana mode localized on the surface the sample at r = L = 20^. Due to tunneling to 
the Majorana core state it also has a negligible energy. 

Fig. 3. 

Excitation energy (in units of the bulk gap A) of the bound states as function of angular 
momentum for several values of the inclusion radius i? at 7 = 0.03. For unpinned vortices 

R = (the bl ue p oints) the spectrum at small / is described well by the semiclassical Volovik 



formula (see 11[), while for higher momenta it gradually approaches the threshold along a 
universal curve that is independent of the inclusion radius R. When the inclusion is present 
the dependence on / is no longer monotonic. It first rises to a maximum at / = — 1 and 
subsequently has a local minimum with energy about that of the I = state. 
Fig.4. 

Dimensionless heat conductance of a single vortex given in Eq.f lTT]) as function 



of inverse temperature in units of A~^. The black line is for unpinned vortex, while the 
corresponding crossover temperature from thermal insulator to conductor for pinned vortices 
on insulating inclusions of radius R = 0.1^ (brown), R = 0.2{ (red), R = 0.4{ (green) are 
much higher. Dashed likes are the contributions of phonons per vortex for B = 0.5T (green), 
0.2T (black), O.IT (blue), length L, = 50nm, density ua = lO^^cm-^, hwo = ^OOK. 
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